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We present measurements of the Nusselt number A/" as a function of the Rayleigh number R in 
cylindrical cells with aspect ratios 0.5 < T = D/d < 12.8 {D is the diameter and d the height). 
We used acetone with a Prandtl number a — 4.0 for 10^ i$ ^ i$ 4 x 10^°. A fit of a powerlaw 
TV — NqW^^^ over limited ranges of R yielded values of 7e// from 0.275 near R = 10'^ to 0.300 
near R = 10^". The data are inconsistent with a single powerlaw for Af{R). For R > 10^ they are 
consistent with TV = aa^^''^^R^^'^ + ha"^^'' R?^'' as proposed by Grossmann and Lohse for a > 2. 
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Since the pioneering measurements by Libchaber and 
co-workers [^2] of heat transport by turbulent gaseous 
helium heated from below, there has been a revival of in- 
terest in the nature of turbulent convection. ^ In addi- 
tion to the local properties of the flow, one of the central 
issues has been the global heat transport of the system, 
as expressed by the Nusselt number TV — Xejf/X. Here 
-^e// = qd/AT is the effective thermal conductivity of 
the convecting fluid {q is the heat-current density, d the 
height of the sample, and AT the imposed temperature 
difference), and A is the conductivity of the quiescent 
fluid. Usually a simple powerlaw 

M = MoW (1) 

was an adequate representation of the experimental data. 

Here R — agd^ IST / kv is the Rayleigh number, a the 
thermal expansion coefflcient, g the gravitational accel- 
eration, K the thermal diffusivity, and v the kinematic 
viscosity. Various data sets yielded exponent values 7 
from 0.28 to 0.31. |^,|| Most recently, measurements over 
the unprecedented range 10^ ^ i? ^ 10^^ were made by 
Niemela et at, and a flt to them of Eq. |] gave 7 — 0.309; 

but even these data did not reveal any deviation from 
the functional form of Eq. ^. Competing theoretical mod- 
els involving quite different physical assumptions made 
predictions of powerlaw behavior with exponents 7 in 
the same narrow range. |^,||-||] Here we mention just 
two of them. A boundary-layer scaling-theory [^,^ which 
yielded 7 = 2/7 2± 0.2857 was an early favorite, at least 
for the experimentally accessible range R ^ 10^^. It was 
generally consistent with most of the available experi- 
mental results. More recently, a competing model based 
on the decomposition of the kinetic and the thermal dis- 
sipation into boundary-layer and bulk contributions was 
presented by Grossmann and Lohse (GL) Q and predicts 
that the data measure an average exponent 7 associated 
with a crossover from 7 = 1/4 at small i? to a slightly 
larger 7 at much larger R. In the experimental range 
the effective exponent -feff = d{ln{N)) / d{ln{R)) , which 
should be compared with the experimentally determined 



7, is only very weakly dependent upon R and has val- 
ues fairly close to 2/7. Thus, it has not been possible 
before to distinguish between the competing theories on 
the basis of the experimental data. 

Here we present new measurements of TV(i?) over 
the range 10^ ;$ i? ^ 4 x 10^° for a Prandtl number 
a = V / K, = 4.0. Our data are of exceptionally high preci- 
sion and accuracy. They are incompatible with the single 
powerlaw Eq. |l], and yield values of 7ey f which vary from 
0.277 near R = lO"^ to 0.300 near R ^ 10^". In particu- 
lar, the results rule out the prediction [^js) 7 = 2/7. For 
R > 10^ a good fit to our results can be obtained with 
the crossover function 

AA = aa-i/i2^i/4^ ^^-1/7^3/7 

proposed by GL |^] for cr > 2. 

We used two apparatus. One was described previ- 
ously 1^. It could accomodate cells with a height up 
to d = 3 cm. The other was similar, except that its three 
concentric sections were lengthened by 20 cm to allow 
measurements with cells as long a 23 cm. In both, the 
cell top was a sapphire disk of diameter 10 cm. A high- 
density polyethylene sidewall of circular cross section and 
with diameter D close to 8.8 cm was sealed to the top 
and bottom by ethylene-propylene 0-rings. Four walls, 
with lengths ranging from 0.70 to 17.4 cm, were used and 
yielded aspect ratios F = D/d = 12.8,3.0,1.0, and 0.5. 
The bottom plate had a mirror finish. It contained two 
thermistors, located 2.5 cm from the center and at 180° 
relative to each other. The fluid was acetone. Its phys- 
ical properties are very well known. | |lO| Measurements 
jm of the convective onset in thin cells yielded values 
of Rc very close to 1708, thus confirming the reliabil- 
ity of the relevant fluid properties. Deviations from the 
Boussinesq approximation, which is usually assumed in 
theoretical work, are relatively small. From Eq. 8 of Ref. 
Ill we estimated a; ~ 1 -h 0.0012AT for the ratio of the 
temperature drops across the top and bottom boundary 
layers. On the other hand, the parameter Q defined by 
Busse iQ, which is usually considered near the onset of 
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convection, can be approximated by Q ~ — 0.05AT and 
does reach significant valued when AT exceeds, say 10°C. 
Our temperature stability and resolution was 0.001°C or 
better. We measured and corrected for the conductances 
of the empty cells, and applied corrections for resistance 
in series with the fluid. The bath temperature was con- 
trolled in a feedback loop and usually fluctuated by no 
more than 0.001°C rms. The bottom-plate temperature 
was also held constant. Its fluctuations were typically de- 
termined by the fluid turbulence, but remained less than 
0.1% of AT. The temperatures of the two thermistors in 
the bottom plate could differ by up to 0.5% of AT due 
to the large-scale flow. The parameters Af and R were 
calculated using the average of the two measurements. 
When only one was used, exponents derived from the 
data typically differed by no more than 4 x 10"''. Usu- 
ally AT was stepped in equal increments on a logarithmic 
scale, holding the mean temperature at 32.00°C. 
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FIG. 1. The Nusselt number as a function of the Rayleigh 
number. Solid squares: F = 12.8; open circles; F = 3.0; solid 
circles: F — 1.0; open squares: F — 0.50: plusses: Ref. | ^ ; 
open triangles: Ref. |l^; solid line: Ref. | |l7[ ]. 



of Liu and Ecke ||T^ for cr = 4, F ~ 1, and a cell with 
a square cross section. The agreement with our data is 
excellent, considering the difference in geometry. 
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FIG. 2. High-resolution plot of the Nusselt number as a 
function of the Rayleigh number. The symbols are as in Fig. m. 
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FIG. 3. The effective exponent as a function of R. Solid 
squares: F — 12.8. Open circles: F — 3.0. Solid circles: 
F = 1.0. Open squares: F = 0.5. Solid line: logarithmic 
derivative of Eq. | with a = 0.326 and b = 2.36 x 10"^. Dot- 
ted line: 7 = 2/7. 



Results of our measurements in four cells of different 
r are shown in Fig. |l| For each cell they cover about 
two decades of R, and collectively they span the range 
of R from 10^ to 4 x 10^°. Of interest here is the range 
R ^ 10^, where one might expect Eq. |^ to be relevant. 
There is a dependence oi Af{R) upon F, as already noted 
by others ||l^. For comparison, we also show in Fig. |l| 
the recent results of Chavanne et al. jl^ for a ~ 0.8 
and F = 0.5 (plusses) and those of Ashkenazi and Stein- 
berg ||l^ for a — I and a cell of square cross section and 
r = 0.72 (open triangles). There is good agreement with 
the former, considering the difference in a. The latter 
are about a factor of 1.8 larger than our results; this dif- 
ference seems too large to be attributed to the difference 
in (T or the geometry and remains unexplained. The solid 
line just above the open circles in the flgure (more easily 
seen in Fig. ||) corresponds to the fit of Eq. |l| to the data 



Figure |l| does not have enough resolution to reveal de- 
tails about the data. Thus we use the early prediction 
7 = 2/7 as a reference, and show /o(7io(A/'i?~^^^) as a 
hmction of logiQ{R) in Fig. ||. If the theory were correct, 
MR^'^/'' should be equal to Aq, i.e. independent of R, 
over its range of applicability. If Eq. |l| is the right func- 
tional form but 7 differs from 2/7, then the data should 
fall on straight lines with slopes equal to 7 — 2/7. Our 
F — 12.8 data (solid squares) are at relatively small R 
and one might not expect Eq. ^ to become applicable 
until R is larger. The F = 3.0 data actually show slight 
curvature, but in any case would yield 7 < 2/7. The 
smaller-F data are clearly curved, showing that Eq. ^ is 
not applicable with any value of 7. In order to make this 
conclusion more quantitative, we show in Fig. ^ effective 
local exponents 7e/ / derived by fitting Eq. ^ to the data 
over various restricted ranges, each covering about half a 
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decade of R. The fits yield values of 7e//(i?) which have 
a minimum near R = 10^. Interestingly, Jeff is within 
our resolution independent of T. 

Next we compare the predictions of GL Q with our 
data. These authors defined various scaling regimes in 
the R — a plane. For a ^ 2, they expect that crossover 
between their regions /„ and ///„ should be observed. 
For that case, Eq. ^ is predicted to apply. One way to 
test this H is to plot y = A/'/(i?^/'*(7~^/^^) as a function 
of X = JJ^/^cr"^/^^. If the prediction is correct, the data 
should fall on a straight line y = a + bx with a and b equal 
to the coefficients in Eq. ||. Our F = 1.0 data are shown 
in this parameterization as solid circles in Fig. ^. The 
solid line is a least-squares fit. One sees that the data 
are fitted extremely well. The coefficients are a = 0.326 
and b = 2.36 x 10~^, in good agreement with the coeffi- 
cients estimated by GL on the basis of other experimental 
data. 
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FIG. 4. Plot of A/'/(i?^/''a~^/") (solid circles) and 
(A/" — 1) /{R^^'''a~^^^^) (open circles) as a function of 
^5/28^-5/84 The straight lines are least-squares fits. IfEq.|| 
is correct, the points should scatter randomly about them. 

We note that in the analysis it was assumed that Af 
is the appropriate variable to compare with Eqs. |l| and 
^. However, one might argue that only the convective 
contribution TV — 1 should be considered. It turns out 
that replacing Af with TV — 1 leads to the same conclu- 
sions. The open circles in Fig. ^ are obtained from our 
F = 1 data when TV — 1 is used. The dashed straight line 
is an excellent fit to the data and yields a — 0.311 and 
6 = 2.58 X 10-3. 

The logarithmic derivative Jeff of Eq. ^ based on 
the analysis of TV (rather than on TV — 1) is shown as 
a solid line in Fig. ^. For R > 10^ it agrees quite 
well with the values determined by local fits of Eq. |l| 
to the data. It would be difficult to demonstrate that 
the small, seemingly systematic, deviations which do ex- 
ist are outside of possible systematic errors in the mea- 
surements. For smaller R, one does not expect the the- 
ory to be applicable. GL estimate that their region /„ 
has a lower boundary below which the Reynolds num- 



ber Re 0.039i?i/V-5/6 < 50. This occurs when 
i? ~ 1.6 X 10^. 

Local fits of a powerlaw to data for TV — 1 are shown in 
Fig. |. There we also show, as a dashed line, jeff based 
on the dashed line in Fig. ^. We see that deviations from 
the prediction Eq. |^ already occur near R — 10^, which 
is an order of magnitude larger than was found in the 
analysis using TV. However, for R ^ 10® the fit of the 
theory to the data is equally good, regardless of whether 
TV or TV — 1 is used. 
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FIG. 5. The effective exponents as a function of R as in 
Fig. ^, but based on using TV — 1 in the analysis rather than 
TV. The symbols are as in Fig. ^. The dashed line is the loga- 
rithmic derivative of Eq.||with a = 0.311 and b = 2.58 x 10"^, 
which corresponds to the dashed line in Fig. H. 
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FIG. 6. Plot of TV/(i?^/''(j-^/i2) ^ ^ function of 
If Eq. H is correct, the points should fall on the 
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straight line, a least-squares fit. The data deviate systemati- 
cally, showing that Eq. |^ is not the right functional form. 

For large a, the theory predicts crossover from to 
IIIu and Eq. |^. For somewhat smaller cr ~ 1, how- 
ever, the crossover should be from to IVu, and the 
prediction then reads [Q 

AA = afT-i/i^i?!/^ + 5i?i/3 (3) 

As a function of a, GL estimate that the transition 
from Eq. || to Eq. ^ occurs near a = 2, but there is 
some uncertainty in this value. Thus, in Fig. ^ we com- 
pared Eq. 1^ with our results for F = 1 by plotting 
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Z/=AA/(i?i/V-i/i2) as a function of :r = R^'^^a^/^'^. 
If Eq. ^ is applicable, this should yield a straight line. As 
can be seen, the data deviate systematically from the fit. 
Thus a transition from /„ to IVu at cr = 4 is inconsistent 
with our data. 
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FIG. 7. Relative deviations 5M/N from fits of (a): Eq. |l| 
(b): Eq. |, and (c): Eq. | to the P = 1.0 data. 

In Fig. ^ we show the deviations 6M /M from fits of 
Eqs. |l| to Eq. | to the T = 1.0 data. For Eq. [l] they 
are systematic, as already expected. For Eq. |^ the fit is 
nearly perfect. We conclude that our data are consistent 
with Eq. |. For Eq. | the deviations arc nearly as large 
as those for Eq. |^. Thus Eq. || is not applicable to the 
data. A similar analysis of the F — 0.5 data over the 
range 10^ ;$ i? < 4 x 10^" gives similar results. 

So far our work has concentrated on measurements of 
Af{R) for acetone with cr = 4.0. We concluded that there 
is no significant range of R over which the powerlaw Eq. |l| 
is applicable, and that the crossover function Eq. ^ pro- 
posed by Grossmann and Lohse provides a good fit to 
the data for R > lO*" where the Reynolds number i?e of 
the large-scale flow is expected to exceed about 50. 
Obviously a great deal of additional high-precision work 
remains to be done, and some of it is indeed under way. 
To provide further tests of the GL predictions, we started 
a systematic study of the Nusselt number A/'(i?, cr) as a 
function of cr, but this work encountered considerable ob- 
stacles. For acetone cr does not vary enough with T, and 
systematic errors in the properties of other suitable 
fluids such as the alcohols prevent a comparison before 
better measurements of A and a/KP are made. Since 
Re is central to the GL theory, its determination is also 
an obvious area for further work. And finally, measure- 



ments of comparable accuracy and precision should be 
extended to the compressed gases where cr ~ 0.7, where 
more direct comparison with much previous work |^,^ is 
possible, and where a different GL crossover function 0| 
should pertain. 
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